Bistability in Feshbach Resonance 
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A coupled atom- molecule condensate with an intraspecies Feshbach resonance is employed to 
explore matter wave bistability both in the presence and in the absence of a unidirectional optical 
ring cavity. In particular, a set of conditions are derived that allow the threshold for bistability, due 
both to two-body s-wave scatterings and to cavity-mediated two-body interactions, to be determined 
analytically. The latter bistability is found to support, not only transitions between a mixed (atom- 
molecule) state and a pure molecular state as in the former bistability, but also transitions between 
two distinct mixed states. 
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I. INTRODUCTION 

The subject of optical bistability [l|, to which Lorenzo 
Narducci contributed greatly during his prime years of 
life, was brought to spotlight again by recent experi- 
mental demonstration of optical bistability in a micro- 
cavity with a cavity field at the level of a single photon 
0, H[. Instead of thermal gases, as typically employed 
by Lorenzo's generation, where the thermal de Broglie's 
wavelength of the particles is far smaller than the inter- 
particle spacing, more recent trend in cavity quantum 
electrodynamics (QED) focuses on cavity systems with 
ultracold quantum gases 0, H[ , opening the door to stud- 
ies such as cavity-mediated bistable Mott-insulator to su- 
perfluid phase transition [H, @ , where the particle statis- 
tical nature becomes an essential feature of the cavity 
problem. In contrast to 1980s, when the surge of in- 
terest in optical bistability was inspired by the prospect 
of its use as a switch in an all-optical computer [lj, the 
current surge of interest in the same topic was, however, 
motivated largely by the equivalence of the cavity con- 
densate system [3| to the cavity opto-mechanical system 
Q. The study of the latter falls into the realm of cav- 
ity optomcchanics, a rapidly emerging field which aims 
to use the cavity-assisted radiation force Q to cool me- 
chanical device, ranging from nano- or micro-mechanical 
cantilevers to macroscopic mirrors in LIGO project, down 
ultimately to their quantum mechanical ground state Q . 

In this paper, we focus on a cavity system containing a 
coupled homonuclear atom-molecule condensate as illus- 
trated in Fig. [T] In the absence of the cavity, molecular 
state 1 2) is coupled only to free atomic state |1) by an 
intraspecies Feshbach resonance Q characterized with a 
strength a' and a detuning e' [the shaded part in Fig. 
QJa)]. In the cavity QED setting, state \2) is coupled, 
besides to state |1), also to excited molecular state |3) 
by a unidirectional ring cavity of a total length C The 
cavity is driven by an external laser of wavenumber k, 
polarization i, and frequency u, which is tuned far away 
from the |3) •<->• \2) molecular transition. In addition, the 
cavity is assumed to possess a sufficiently large intermode 
spacing C/c (with c being the light speed in vacuum) so 
that only the cavity mode with a longitudinal mode fre- 
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FIG. 1: (Color online) (a) The schematic energy diagram 
(with the shaded part representing the intraspecies Fesh- 
bach process) for (b) a unidirectional ring cavity + coupled 
homonuclear atom-molecule condensate system. 



quency u c closest to oj is relevant to our study, where 
uj c £/c equals integer multiple of 27r. Further, the cav- 
ity mode is assumed to overlap with the condensate in a 
spatial region, characterized with a length L and a cross- 
sectional area A, large enough so that the condensate can 
be treated as a uniform system with an effective volume 
V a = LA and a total atom number N a (counting those 
in molecules). 

A similar cavity + condensate system has recently been 
studied (Tol j , in which atoms are converted into molecules 
by photoassociation as opposed to magnetoassociation 
(Feshbach resonance) in our model. There, absorption 
of a cavity photon will convert two atoms into a molecule, 
and emission of a cavity photon will dissociate two atoms 
into a molecule. This is reminiscent of the absorptive 
bistability in a driven optical cavity containing an en- 
semble of two-level atoms. In contrast, the cavity field in 
our model is to introduce a phase shift to the molecular 
field, not the population exchange between atoms and 
molecules. As a result, the bistability in our model is 
of dispersive nature. Further, interaction between elec- 
tronic dipoles and a cavity field in Ref. [lfj gives rise to 
atom-molecule coupling, which, since the cavity field is 
itself a dynamical variable, changes with time. In con- 
trast, in our model, it is the hyperfine interaction - the 
coupling between electron spins and nuclear spins of two 
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colliding atoms that result in atom-molecule coupling, 
which is therefore fixed by the atomic internal structure, 
independent of cavity field. In this respect, our model is 
analogous more closely to a cavity s yst em with a spin-1 
condensate (also recently propose d II 11 1 than to a cav- 
ity system with photoassociation [lOj; atom exchanges 
among different internal states in the spin-1 condensate 
are accomplished by spin-exchange interaction, which is 
also independent of cavity field. 

The interest in Feshbach resonance stems primarily 
from its use as an effective tool to coherently create 
molecular BECs from the existing atomic BECs. In 
a typical experiment, conversion of atomic BECs into 
molecular BECs is carried out by ramping the Feshbach 
detuning across the resonance. This method relies on 
the existence of a mixed (or dressed) atom-molecule state 
[lH, and the ability of this state to change its composi- 
tion from predominantly atomic to predominantly molec- 
ular species when the Feshbach detuning e' is tuned from 
above to bellow the resonance. The question that we 
want to pursue, in this paper, is how the molecular pop- 
ulation in state |2) can be made to vary with the Feshbach 
detuning in a bistable fashion, instead of monotonously 
as in typical situations. A bistable crossover adds new 
meaning to the Feshbach resonance: whether the system 
is in the atomic or molecular extreme is determined not 
only by the Feshbach detuning but also by the history of 
the system. 

To this end, we first formulate, in Sec. II, a semi- 
classical mean-field description of our model in the limit 
where both the excited molecular field and the optical 
cavity field can be adiabatically eliminated. We then 
explore the matter-wave bistability due to the two-body 
s-wave collisions in Sec. Ill A and that due to the 
cavity-mediated effective interaction between two Fes- 
hbach molecules in Sec. Ill B. Finally, we provide a 
summary in Sec. IV. 



II. THE BASIC EQUATIONS 

In this section, we take a semiclassical approach, in 
which optical fields are treated classically while matter 
fields are treated quantum mechanically, to formulate a 
theoretical description of the proposed cavity + conden- 
sate system. This is the same approach that Lorenzo 
embraced in many of his works [13 1, except that quan- 
tization is now performed at a level for a many-body 
system, instead of a single-body system as in a typical 
semiclassical laser theory. To begin with, we expand op- 
tical field E (r, t) in terms of the slowly varying amplitude 
F in space r and time t according to 



E (r, t) = ^Fee l{kz - 



and matter field ip ( r i t) in terms of ipi — Ci/V N a accord- 
ing to 



i, 1 \l) + ^ 2 \2)+^ kz -^\^)\, (2) 



where Cj denotes the operator in momentum space for an- 
nihilating a bosonic particle in state \i) and n a — N a /V a 
is the total atom number density. The expansion in 
Eq. @ is carried out in a frame rotating at the laser 
frequency u>. In arriving at Eq. ([2]), we have assumed 
that the particles in states |1) and |2) are all condensed 
to their respective zero momentum modes, and those in 
state 1 3) to the mode with hk momentum in accordance 
with momentum conservation during photon-atom inter- 
action. 

The coupled atom-molecule system, within these ap- 
proximations, is then described by the following Hamil- 
tonian 



H/hN a = e'i)\i) 2 + I y V^i + h.c 
+ n a ^^\i>]^i)i = 1 or 2) 



?,U,2 



+ (e' - A Q ) flfa - (fi$$, 



+ h.c), (3) 



where repeated indices are to be summed from 1 to 2. In 
Eq. (j3|), the first line describes the Feshbach resonance 
of strength a' , the second line the s-wave collisions of 
strength x'ij (= Xji) between states \i) and \j), and the 
last line the part of Hamiltonian involving excited state 
|3). In the last line, the first term denotes the energy of 
state 1 3) in the rotating frame, where A a is the laser de- 
tuning, and the second term stands for the laser-induced 
electric dipole interaction, where il = [i-^F jlh is the 
Rabi frequency, /132 = (3|/t • e\2) the matrix element, 
and \x the electric dipole moment operator. Finally, col- 
lisions involving the final excited state |3) are ignored 
since state |3) remains virtually empty in our model. 

An important concept in the semiclassical approach 
is the macroscopic polarization defined as P (r, t) = 

^V>t (r,t) (lip (r,i)\. This polarization is found, when the 

use of both Eq. ((2]) and the selection rule are made, to 
possess the same mathematical form as the optical field 
in Eq. ©, 



P(r,t) 
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(4) 



where P = 2n a ^2 ^V^V'I ) represents the slowly varying 
part of the polarization. The evolution of the optical 
field is then governed by the Maxwell's equation, which, 
under the slowly varying envelope approximation, takes 
the form 



dz 



on 

~dt 



(5) 



c.c, 



(1) 



where [1q (eo) is the magnetic (electric) permeability in 
vacuum. Equation ([3]) clearly shows that polarization 
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plays the role of a bridge between the classical optical 
field in Eq. (JTJ) and the quantum matter fields in Eq. 

©■ 

The final component in the semiclassical approach per- 
taining to any QED problems is the boundary condition, 
which, in our case and under the assumption that both 
input and output mirrors have the same reflectivity R 
(transmissivity T = 1 — R), can be cast into the form 

m 



represents the phase shift to the molecular field induced 
by the cavity optical field, where N c — Y 2 /g 2 repre- 
sents the number of injected intracavity photons and 
rj = g 2 / A a measures the effective atom-cavity coupling 
strength. (In this paper, without loss of generality, we 
limit A a > so 77 is always positive.) 



III. DISCUSSIONS 



fi (0, t)=TY + RQ (L, t - At) e 



■;a c £/c 



(6) 



where Y is the scaled amplitude of the incident field, 
At = (C — L) I c the transit free propagation time inside 
the cavity, and A c = w — ui c the cavity mode frequency 
detuning. Equation ((6]) links the field entering {z = 0) 
to that leaving the condensate [z = L), and hence imple- 
ments the concept of feedback, which is the most impor- 
tant feature of any optical cavities. The school led by 
Bonifacio, Lugiato, and Narducci distinguishes itself by 
its insistence on a rigorous (first-principle) treatment of 
the boundary condition - a treatment made up of both 
a transformation mapping two non-isochronous events in 
Eq. ([5]) into two isochronous events and a set of 
conditions embodying the notion of "mean-field limit" 
16]. Following this treatment, we explicitly eliminate 
the spatial derivative in Eq. ([5]), rendering Eq. ([5]) into 



— = (iA c - K)ll- 
at 



kY 



W 2 N a 



(7) 



In our system, stationary solutions can be divided into 
mixed atom-molecule states in which each species has a 
finite density nj, and pure molecular or atomic phase in 
which one of species has a zero density. In this section, 
we study bistability in the Feshbach process by focusing 
on the mixed atom-molecule state where each species can 
be described by the field = ■ v /n7e l6 ' i characterized with 
a well defined phase 9i. (Note that such a description 
cannot be applied to pure molecular or atomic phases as 
the phase associated with an empty component is not 
defined.) For such a state, we can apply the total atom 
number conservation n% + 2n2 = 1 and simplify Eqs. ([5]) 
into a set of equations involving only two variables - a 
molecular density n = 712 and a phase mismatch 8 = 
82 — 28i. This set of equations, in a unit system in 
which e = (e' — X12 — 2xn)/a is defined as the effective 
Feshbach detuning, x = [X22 + 4 (xn ~ Xi2)]/a as the 
effective Kerr nonlinear coefficient for the molecular field, 
S = A c /k as the cavity detuning, and r = at as the time, 
take the form 



where n — c\hiR\/C is the cavity damping rate, g = 
H32^/oj/2he V c the "Rabi frequency" per photon, and 
V c = LA the total cavity volume. 

In this paper, we restrict our study to the parameter 
regime in which both the time scale for the optical field 
(on the order of 1/k) and that for the molecular field 
of state |3) (on the order of 1/A Q ) are far shorter than 
those for the Feshbach degrees of freedom, namely, the 
matter fields corresponding to states |1) and |2). This 
restriction allows us to adiabatically eliminate the former 
fast variables in favor of the latter slow ones, reducing, 
from the Heisenberg's equations for ipi and the Maxwell's 
equation for ft, a set of equations involving only the slow 
degrees of freedom ipi and 1^2 '■ 

i ~aT = XliUi ^ 1 + a ^ 2 ^i' ( 8a ) 
i d ^ = [e' + X2 l n l + F(n 2 )]^2 + ^l. (8b) 

where we have adopted the standard mean-field theory, 
treating ^ as c-numbers fy. In Eqs. (O, n\ — |V , i| 2 I 
n-2 = 1^21 , a = a'y/n^, and Xij = x'ij n a are the renor- 
malized quantities, and 



F(n 2 ) = 



K 2 T]N C 



k 2 + (A c - 7]N a n2) 



(9) 



dn 
dd 

— = e + X n+ -- 
or 2 x/n 



2n) y/n sin 8 
1 1 - 6n 



cos 8 + f (n) , 



where 



f(n) 



B/C 



1 + (S - CnY 



(10a) 
(10b) 

(11) 



In arriving at Eqs. (flU)). we have replaced N a and N c in 
Eq. (|9]) in favor of two unitless parameters, 



and 



C = T)N a / k = g 2 N a /A a n, 



B = V 2 N c N a / (2kq) 



(12) 



(13) 



In contrast to C, which is bose enhanced by atom num- 
ber only, B is bose enhanced by both photon and atom 
numbers. Note that when A Q is replaced with the decay 
rate of the excited state |3), C becomes the so-called 



atomic "cooperative parameter" [17J. As we will see 
shortly, cavity-mediated bistability depends crucially on 
the values of C and B. 

As in other studies [H, [H], Eqs. ([TDJ of the type 
including Feshbach resonance support two branches of 
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steady-states: one with 9 — and the other with 8 = it. 
This feature is expected since at zero temperature the in- 
traspecies Feshbach resonance represents a matter wave 
analog of the second harmonic generation in nonlinear 
optics where the phase-matching condition plays an im- 
portant role. In this paper, without the loss of generality, 
we take a > 0. Under such a circumstance, the branch 
with 9 = 7r not only always has a lower energy than the 
branch with 9 = 0, but also has the property of consist- 
ing primarily atom species in the limit of a large positive 
Feshbach detuning. For these reasons, we will focus on 
the branch with 9 = n, determined at steady state by 



1 1 - 6n , , 

-Xn+ - — ■= / (n 

2 y/n 



(14) 



where for notational simplicity, same symbols are used 
to stand for the steady state variables. 

To carry out the stability analysis, we apply the stan- 
dard linearization procedure that Lorenzo used exten- 
sively in the context of his interest in laser instabilities 
[20l l2~i| . and derive from Eqs. (|10p a set of linearly cou- 
pled equations 



d_ 



6n 



(1 - 2n) ^ 

an 



5n 
68 



(15) 



where 5n and 89 are small departures from the corre- 
sponding steady state variables, and e is given by Eq. 
(fT4]) . The eigenvalues of Eqs. (|15l) are then found to 
take two values: iy — (1 — 2n) ■y/nde/dn, which, since 
n < 0.5, clearly indicate that only when de/dn > or 
equivalently when dn/de > can the eigenvalues become 
complex. In another words, any state at which the slope 
of n as a function of e is positive is unstable. Note that 
such a conclusion may not hold, if we do not impose, in 
the previous section, the conditions that allow ^3 and fl 
to be adiabatically eliminated. In the case of absorptive 
optical bistability, it is well known that the upper branch 
which is stable according to the slope criterion may be- 
come unstable; instability there is manifested in the form 
of self-pulsings [22|. 

Thus, we see that the stability analysis here amounts 
to analyzing the points at which de/dn = 0, which, by 
definition, simply corresponds to the critical transition 
points in a bistable (or multi-stable) system. As a 
result, in what follows, we carry out bistability study by 
focusing on the equation 



d(n)=x-h(n), 



(16) 



derived from Eq. (|14[) under the condition that de/dn 
0, where 



d(n) = —AB 



(5 - Cn) 



h(n) 



I + {6- Cn)' 



i n -l/2 _ 1 -3/ 2< 



(17) 
(18) 



A. Collision-Induced Bistability 

In our model, the molecular field (at state |2)) is sub- 
ject to two types of self phase modulations: one, de- 
scribed by a Kerr type of nonlinear term x n > originates 
from short-range s-wave scatterings, and another, de- 
scribed by / (n) in Eq. (Ill[) . stems from cavity- mediated 
long-range two-body collisions. In order to differentiate 
their roles in the formation of matter wave bistability, we 
first remove the cavity component and study the bista- 
bility due solely to the Kerr nonlinearity by solving 



X = h(n), 



(19) 



obtained from Eq. (flU]) by setting d(n) = 0. As one 
may easily verify, h (n) in Eq. (|18p is a monotonously 
increasing function of n. As a result, in order for Eq. (fT9|) 
to have real roots within n < 0.5, x must be less than 
h(n = 0.5) or 



X < Xth 



-2V2. 



(20) 



This condition is similar to that of Ref. [2j| for a het- 
cronuclear atom-molecule system with an interspecies 
Feshbach resonance [l8(. It simply reflects the fact that 
for bistability to take place, there must be a sufficiently 
strong positive feedback between the molecular popula- 
tion and the effective Feshbach detuning e + x n - A nega- 
tive x fulfills this positive feedback; as can be seen, with a 
negative x, an increase in the molecular density decreases 
the effective detuning, which, in turn, further increases 
the molecular density, or vice versa. Such a chain of 
positive reaction under condition ([20)) can lead to the 
formation of the critical transition points around which 
the molecular population changes in a runaway fashion. 
Indeed, under condition (|20l) . Eq. (fl6|) is found to sup- 
port a critical point with a critical molecular density 
given by, 



1 



16 X 



X 



16 x 



X 



1 

2x 



-^r, (21) 



and a critical Feshbach detuning e cri determined by Eq. 
(IT41 when n is replaced with n^ r { in Eq. (|2"Tj) . Figure 2(a) 
shows a typical example of bistability (with x = %Xth)- 
In addition to the mixed state, there is a pure molecular 
state with n = 0.5 obtained from Eqs. (jHJ), not from Eqs. 
(|10p which, by definition, only holds for mixed states. 
(Note that the pure atomic state is prohibited by the 
nature of intraspecies Feshbach resonance [lj| H3, HH.) 
The interception between this pure state and the mixed 
atom-molecular state defines the second critical point 



,(2) 



5 f (2) 



X 
2 



y/2. 



(22) 



5 






0>) 






A 




en \ 





0.5 



1.5 




0.2 0.3 0.4 0.5 



0.5 = 
0.4 
0.3 

0.2^- 
0.1 

— 



(b) 



(-0.58, 0.2) 



FIG. 2: (a) Molecular population n as a function of e when 
X = 2%th = —4y/2. (b) e^' and at different values of x- 



The threshold for bistability is reached when the two crit- 
ical points become degenerate. Clearly, this happens at 
n = 0.5 when x = Xth- Figure HJb) shows that when x 
is bellow Xthi the size of the hysteresis loop (measured 



e^.;-) increases with |x|- 



by £ (2) - e (1) . 

To see the implication of condition (|20[) to a realistic 
system, consider, for example, the Feshbach resonance 
located at the magnetic field 85.3 mT in 23 Na (2ol |27|. 
This resonance has a width of 0.95/iT or equivalently a 



Feshbach coupling strength of a' = 4.22 x 10 



-3/2, 



S . 

Taking the total atom number density to be 10 20 m~ 3 
and using 3.4 nm as the s-wave scattering length for 
sodium atoms [28|, we find that xn = 1.18xl0 4 s _1 
and a =4.22 xlO 4 s _1 = 3.58xii- If we further as- 
sume X22 = Xil] we see from Eq. (|20[) that this requires 
Xia > 3 (X22 + 4xii + 2V2a) > 3.77xil 



B. Cavity-Mediated Bistability 

The above example means to illustrate that Eq. (|20|) 
can be fulfilled only when both the total atom density 
and the two-body interspecies collisional strength are suf- 
ficiently large, a condition which is difficult to meet under 
typical systems. In this subsection, we turn our atten- 
tion to the cavity model in Fig. [TJ and pursue, from Eq. 
(1161) . the question of under what cavity parameters can 
bistability occur even when condition (1201) breaks down. 
Note that unlike Eq. (TT51) . which only contains one crit- 
ical point, Eq. (|16[) can typically support two critical 
points. Thus, the threshold for bistability in this cavity 
model can, in principle, occur at any value of n, instead 
of always at n = 0.5 as in the bare Feshbach model dis- 
cussed in the previous subsection. To proceed, we first 
change Eq. (|16[) into a quartic equation for^/y 



(Vv) 4 -7^+1 = 0, 



(23) 



where y — Cn — 5 and 7 = 2y/B/ [x — h (n)}. This 
change of variable is motivated by the realization that 
when condition (|20[) breaks down or equivalently x + 
2-\/2 > 0, the quantity \ — h (n) is always positive so 



FIG. 3: (Color online) (a) B t h (n) as a function of n. (b) 
Molecular population n as a function of e for B — B t h (0.2) = 
4.73 (solid line), B = 2B th (0.2) (dashed line), and B = 
4:B t h (0.2) (dot-dashed line) for a cavity system with C = 20, 
X = 0, and 5 = 8 th (0.2) = 3.42. 



that only when Cn — S > can Eq. ([16| holds. To esti- 
mate the threshold condition at a given n, we regard Eq. 
(|23p as a transcendental equation for S, and require Eq. 
(|23| to support a real root of multiplicity 2 (another two 
roots are a complex conjugate pair). This requirement 
allows us to conclude that bistability develops when 

B>B th {n) = ^j=[ X -h{n)\, (24) 

and the bistability threshold at a given n (and C) is 
reached when B = B t h (n) and S = 5th i n ) = Cn — 1/ v3- 

B t h (n) as a function of n is displayed in Fig. EJa). 
This simple threshold relation seems to hold quite well 
as long as C is sufficiently large. Consider, for example, 
a cavity system with C — 20. Figure |3tb) shows how 
n changes with e under different values of B. The solid 
line, produced with B = B t h (0.2) = 4.73 [point A in Fig 
[3ja)] and 5 = 5th (0-2) = 3.42, shows that the bistability 
threshold indeed occurs at the theoretically predicted lo- 
cation with n — 0.2 [and e = —0.58 determined from Eq. 
(TT4")) ]. The dashed and dot-dashed lines in Fig. EJb), 
produced with B = 2B th (0.2) and B = AB th (0.2), re- 
spectively, clearly shows that the usual bistable behavior 
emerges when B is increased beyond its threshold value 
B t h (0-2). Recall that without a cavity, a bistable tran- 
sition can only take place between a mixed and a pure 
molecular state. In the present situation with an opti- 
cal cavity, we see that a new feature appears - a bistable 
transition can also take place between two different mixed 
states (the dashed line). 

To see what condition means to the cavity field, 
we consider a small-sized micro-ring cavity of total length 
C = 200 fim and finesse T = 3.14 x 10 5 . The cav- 
ity is driven by an external laser tuned A a — 2tt x 100 
GHz away from the D2 line, characterized with a wave- 
length 780 nm and a linewidth 27r x 3 MHz, of sodium 
atoms confined to an effective spatial region of L = 30/im 
and A — (10/im) . In such a cavity-condensate sys- 
tem, we have 77 = 2n x 0.0052 MHz and k = 2n x 23.87 
MHz. Then, by using the same Feshbach resonance 
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and the atom number density in the previous subsection, 
we estimate that the minimum number of photons that 
must be present inside the cavity to produce bistability is 
8.6 x 10 3 . This figure is three orders of magnitude smaller 
than the threshold photon number in a typical laser [l7| , 
and can be further reduced with an appropriate choice of 
the system parameters. 

Matter wave bistability at small photon numbers can 
be understood as follows. In contrast to the phase shift 
due to the s-wave scattering, which is linearly propor- 
tional to the molecular density, the phase shift, arising 
from the feedback between optical and matter fields, is 
nonlinearly proportional to the molecular density by way 
ofEq. (fTTj) in a resonant fashion. On one hand, the sensi- 
tivity (the change of this phase shift versus the change of 
the molecular density) is bose-enhanced by the collective 
nature of the condensate system. On the other hand, the 
peak of this phase shift for a given cavity photon number 
can be significantly amplified in a microcavity environ- 
ment where photons are confined into a tiny volume. As 
a result, molecular bistability is possible under a weak 
cavity field. 

IV. CONCLUSION 

In this paper, we have studied the matter wave bista- 
bility in an intraspecies Feshbach resonance model with 



and without the assistance of an optical cavity. In par- 
ticular, we have arrived at a set of conditions that allow 
the bistability thresholds under the two different settings 
to be estimated analytically. In the absence of a cavity, 
bistability is possible only when the effective Kerr nonlin- 
earity stemming from s-wave scatterings [Eq. (|20[) ] is suf- 
ficiently negative. In the presence of a cavity, even when 
condition (f!HI)) breaks down, bistability can still occur, 
provided that B in Eq. (|13l) , a key parameter describing 
the cavity-mediated two-body interaction, is sufficiently 
large [Eq. [M] . An important difference between systems 
with and without a cavity is that the former supports one 
stable mixed state while the latter can support two dif- 
ferent stable mixed states. (Both systems contain a pure 
molecular state.) Thus, a bistable transition in the lat- 
ter system can take place not only between a mixed state 
and a pure molecular state as in the former model, but 
also between two different mixed states. 
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